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The production of r) mesons in NN collisions in the near threshold region is analyzed with special 
emphasis on the role of final state interaction which is calculated using three-body scattering theory. 
The three-body aspects are shown to be very important in interpreting experimental results for the 
total cross section of pp rjpp. Furthermore, energy spectra and angular distributions of the 
emitted particles are calculated and compared to experimental data. 

PACS numbers: 13.60.Le, 21.45.+V, 25.20.Lj 



O ■ I- INTRODUCTION 

o ; 

• The production of 77 mesons in two-nucleon collisions depends strongly on the final state interaction (FSI) as has 
^ ' been stressed in previous work X^t^^^^ Notwithstanding this fact, the interaction between the final particles 

was treated in quite an approximate way. Only the NN interaction has been considered in a more precise manner. 
The interaction of the 77 meson either was not considered at all or simulated by means of a Watson enhancement 
0^ factor. In view of this apparent lack of a correct treatment of FSI we want to study this reaction by using a rigorous 
^ ■ three-body model. 

The properties of the rjNN system, containing three free particles in the asymptotic region, were investigated in 
where FSI effects in r] — d reactions as well as in 77 photoproduction on a deuteron were studied. As was shown in 
^^<-^ [ detail, the three-body aspects of FSI are very important, at least near the threshold region where S waves give by far 

• the largest contribution to the production process. It is obvious that the three-body S wave part is most sensitive to 
' FSI in view of the dominance of the S wave interaction in the NN and r]N subsystems at low energies. Furthermore, 

it is intuitively clear that primarily the central region of the system, where S waves dominate, is distorted by multiple 
scattering events. As a consequence, close to threshold the rjNN interaction strongly enhances the ■fd rjnp cross 
section. At higher energies, where the large spatial extension of the deuteron comes into play, an essential fraction of 
higher partial waves contributes increasingly to the r] production mechanism. This rapidly reduces the relative S wave 
contribution to the reaction amplitude, so that the role of FSI diminishes considerably with increasing energy. As a 
result, already in the region of cm. excess energies above 100 MeV the behavior of the cross section is well explained 
^ ] in the plane wave approximation for the final state Q. 
^ . Turning now to the process NN — > riNN, firstly we would like to note that in this case the influence of FSI must 
^ ■ be even more important. Here the characteristic radius R of the reaction, determined hy R 1/q with q as typical 
• w^ , momentum transfer, is only about 0.3 fm. This value is much smaller than in photoproduction, where not too close 
^ ■ to the threshold the characteristic range is given by the deuteron radius r^, i.e., i? ^ w 2 fm. Therefore, one may 
M , expect that for energies for which pR < 1 is valid, with p denoting the maximal cm. momentum of the final particles, 
- - I the final state will be dominated by S waves, for which the role of the rjNN interaction is crucial. Indeed, as will be 
shown below, in the region from threshold up to an excess energy Q = 60 McV, the contribution of higher partial 
waves in the final state is of much minor importance compared to the S wave contribution. In this connection, the 
reaction NN —^ ijNN seems in general to be more suited for studying the interaction in the r]NN system than 77 
photoproduction, where FSI has a smaller influence. This advantage is however strongly reduced by a poor knowledge 
of the 77 production mechanism in NN collisions. 

In the present paper we investigate the reactions np — > r]np and pp — > rjpp in the region of excess energy up to 
Q = 60 MeV which is covered in recent experiments as well as in theoretical studies (see e.g. and references 
therein). The two questions which we address here are: 

(i) T o what extent are the three-body aspects essential for a quantitative description of FSI? As was already stressed 
in TO] , the importance of a three-body treatment of FSI is dictated by the proximity of a virtual pole in the rjNN 
amplitude to the physical region Tlj. Because of this property a perturbative approach to the rjNN dynamics, where 
only the interactions in the two-body NN and j]N subsystems are taken into account, is of little use. In the present 
study we would like to demonstrate this fact for 77 production in NN collisions. 
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(ii) What is the role of higher NN partial waves? According to Nakayama [l^l , the puzzling second bump in the 
pp invariant mass distribution of pp — > r]pp measured at Q — 41 MeV in could be understood as arising mostly 
from a contribution of the ^Pq state of the final nucleon pair. However, the size of P wave contribution appears to be 
comparable with that of the ^Sq one (see e.g. Fig. 4 in ^^l). Such a large P wave effect in the region where the NN 
cm. kinetic energy is limited to Tjvat < Q = 41 MeV (and in fact the maximum of the distribution of Tj^i^ is located 
at even a considerably lower value because of the preference of higher rj momenta) is quite surprising. Therefore, we 
have investigated again the influence of S, P, and D partial NN waves in both isospin states using the Paris NN 
potential. Another point, closely connected to this, is the role of the short range behaviour of the NN interaction 
which is enhanced by quite large momentum transfers accompanying the 77 production. As will be shown below, this 
leads in particular to an important role of the tensor force in the final np triplet state for np ijnp because the ^Di 
state, excited in the r] production process, couples to the '^Si partial wave during subsequent NN rescattering. 

First we review briefly in the next section the most essential formal ingredients on which our calculation is based. 
For the driving NN — > NN* transition potential we take a meson exchange model including tt, p, and 77 exchanges. 
The two-body rjN and NN interactions are parametrized in terms of separable t matrices. The dominant 5' wave part 
of the rjNN wave function is calculated using three-body scattering theory. In Sect. lIIII our results are compared with 
experimental data and in Sect. Hvl the most important conclusions are summarized. In an appendix we list detailed 
expressions for the NN — > NN* amplitude. 



II. FORMAL INGREDIENTS 



In the center of mass system the reaction to be studied reads in detail 

N{E,p) + N{E, ^p) ^ N(E^,n) + N{E2,P2) + vi^^, q) , (1) 

where the energies and momenta of the participating particles are explicitly defined for later purpose. The reaction 
is governed by the transition matrix element 

Mtai^,S'M'^,sais = {q;Pi,P2;S'Msr,TMT\M\SMs;TMT,p) , (2) 

where initial and final states are characterized by the total spin S (S") and isospin T of the nucleon pair. The 
calculation is done in the region of cm. excess energy Q up to 60 MeV. Thus the use of nonrelativistic expressions 
for the energies of the final particles is justified, whereas for the initial state relativistic expressions are taken. 
The cm. cross section for pp — > rjpp is then given by 

where the relative momentum of the two final nucleons Pr = [pi ^ P2)/2 is introduced. The factor j results from the 
average over initial spin states. The corresponding np — > rjnp cross section is obtained from (PJ by the substitution 

Mii^S'M'g,SMs 2 $Z Mto,S'M'^,SMs ■ (4) 
T=0,1 

For the rj production mechanism we have adopted the usual model where first an N* resonance is excited on one of 
the nucleons which subsequently decays into the rjN channel. The submatrix Mnn" associated with the transition 
NN ~> NN* is expressed as a linear combination of eight pseudoscalar operators (j = 1, 8) built from the vectors 
p, p' , and spin operators ai and (72 

8 

Mnn' {p,p') = X! ^jiP'P'' ^1,^2) , (5) 
i=i 

where p' denotes the relative momentum of the NN* system. The amplitudes Mj are functions of only p = \p\, 
p' = \p'\ and p- p', and their specific form depends on the model for the 77 production mechanism. For the operators 
we have chosen the following forms 

Si(2) = 

^3(4) = CT2 • _p , 
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^5(6) = ^("^1 ^ ■P^'\ 

S7(8) = ^\/5{(al®<T2)''l®(p^'^®(pxp'))^''}™ ■ (6) 



The latter two operators are tensor operators in two-nucleon spin space describing a spin transfer of 2 in the transition 
to the NN* state. For the impulse approximation (lA), where FSI is ignored, the amplitude Mtmt,S' m'^,sMs © is 
given by 



-{-lf+^(AdNN'{vJi) - {-lf+^MNN'{-p,Pi))FN'{q,P2)\ \SMs;TMt) . (7) 



Here the function FN*{q,p) describes the propagation of the NN* configuration followed by the decay N* rjN. 

With respect to the spin structure of the amplitude, only three types of transitions are possible, namely S = 
S' — 1, S = 1 S' — 0, and S" = 1 5" = 1. The transition 5 = — > 5' = is excluded since it does not conserve 
parity. It is easily seen that the operators Si through S4 contribute to all three types of transitions. The operators 
E5 and Sg, being antisymmetric under exchange of ai and (T2, appear only in channels with a change of spin, i.e. 
S ~ ^ S' ^ 1 and 5 = 1 ^ 5" = 0, whereas S7 and Ss, which require flipping of the spins of both nucleons, 
contribute only to 5* = 1 5' = 1. 

At low energies, the channel 5 = 1 — > S" = is strongly favored in the pp rjpp reaction since it allows 
the final nucleons to be left in a relative S wave whose contribution is further increased by FSI. For the reaction 
np — > rjnp both spin changing transitions are possible near threshold. Apart from a trivial isospin factor 3 appearing 
for T — 0, their relative importance depends crucially on the interference between tt and p exchange 2j. The transition 
5 = 1 ^ S" = 1, which requires either the nucleons to be in a relative P state or the 77 meson to be emitted at least 
in a P state is suppressed at low energies. 

After these more formal considerations we will now turn to the outline of the rj production model used in the present 
work. The 77 production mechanism via N* excitation through meson exchange was considered in a variety of previous 
studies QjHSSIS- In the present work, we have adopted this approach including only tt, p, and 77 mesons whose 
contributions was investigated in greater detail than those of w, a and other mesons. The corresponding Lagrangians 
are taken in the customary form 

Ct,nn = ig.„NN FT,[q^) u{p')"f^ u{p) T ■ cj)^ , 

CrrNN' = g^N N' F* {q^ j^) u{p') u{p) T ■ (j)^ , 

CriNN = igriNN Friiq^) u{p')^^ u{p)(l)n , (8) 

C,^NN' = gr,NN' F*{q^j^) u{p) u{p) (f>rj , 

CpNN = Fp{q^)u{p')l^gv7'' + i^^<^'"'q.^u{p) e,f-$p, 

CpNN* = iFp{q^) u{p') j^j^^'^ j^j^ cr^'^qi' u{p ) e„ t ■ $p , 

where q ~ p — p' denotes the four momentum of the exchanged meson and q^N — {^Inq — TT^aP) / {Mj^ + ma) the 
relative meson nucleon three momentum with a € {tt, rj}. Further, each meson nucleon vertex aNN with a G {tt, 77, p} 
contains a form factor of monopole type 

F^iq^)^f^^^. (9) 
The corresponding aNN* transition vertices with a e {tt, t]} are regularized by form factors of the form 

Ki'i!N)-jJ^- (10) 

For the pNN* vertex the same form factor as for pNN is used. With this meson exchange model the functions Mj 
appearing in (jSJ can be written as a sum over the participating mesons 



(11) 
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Detailed expressions for M a £ {tt, t], p} are listed in the appendix. 

Following the work of the main calculation is based on the assumption that the NN NN* transition is 
dominated by p exchange. The contribution of tt meson exchange to the amplitude is about three times smaller. The 
rj meson with an t]NN coupling g^jqj^/'^TT = 0.4 from provides only a tiny correction at the level of about 3 %. 
The pNN* coupling constant was treated as a free parameter. It was adjusted in such a way that the total cross 
section in both channels, pp and np, is normalized to the corresponding data with the same normalization factor 
C„ < 1. This factor takes effectively into account the damping of the initial beam via inelastic collisions other than 
NN r]NN. 

The parameters appearing in JHJ to Hl()|l are listed in Tables HI and Hll Those for the aNN couplings, used also in 3 
are give n in except for g,fNN which was taken from H^f. The ttNN* and rjNN* parameters are taken from [lq| 
and [13 as is explained below. With this set of parameters we obtain as normalization factor C„ = 0.18, which is 
appreciabl y sm aller than the value C„ ~ 0.55 of (3| and also less than the values 0.23 and 0.2 obtained respectively 
in 01 and In view of some obvious limitations adopted here for the production potential, e.g., neglect of direct 
emission of an rj meson by one of the nucleons as well as the contribution of other mesons, we do not put too much 
significance on this quite small value for the damping factor C„. 

Clearly, the present model is only one of possible ways to parametrize the NN NN* transition amplitude. As 
was discussed in , the isovector exchange with a dominance of p over tt exchange is consistent with the observed ratio 
of about 6.5 of the np — > r]np to pp — s- rjpp cross section '18^ as well as with the characteristic (a — cos^ 0,,)-dependence 
of the rj angular distribution ^IQj . For this model the transition with T = strongly dominates over T =^1 and, among 
others, results in a large probability of deuteron formation 2] in agreement with experimental results ^ISj. However, 
as has been shown in ^ and ^ , more sophisticated treatments of the initial state interaction as well as inclusion of 
genuine two-body meson production mechanisms favour tt dominance. Therefore, in order to broaden our theoretical 
basis, we consider in addition the second type of NN — > NN* potential, where the p meson is neglected. Another 
purpose of using such an extreme approximation is to see to which extent the FSI effects depend on the details of 
the NN — > NN* transition. One might anticipate the result that due to the short-range nature of the NN NN* 
interaction such a dependence can not be very significant. However, in view of the basic theoretical uncertainties with 
respect to the production mechanism, this question has to be investigated quantitatively. Indeed, as will be shown 
below, the sensitivity of FSI effects to the variation of the relative strength of tt to p exchanges in the production 
potential, although not very dramatic, results in a variation of FSI effects of about 30 %. 



III. RESULTS AND DISCUSSION 



Our method of solving the three-body equations is based on separable representations of the driving two-body 
forces. For the rjN interaction we use the isobar model of [l^ . where the resonance parametrization of the partial 
wave amplitudes is fitted to the ttN —^ ttN, ttN — > ttttN, and Tr~p rjn channels. In particular in the Sn channel 
only the 5ii(1535) resonance was taken into account. Apart from the original set of Sn parameters presented in 
|16l |. which gives for the r]N scattering length Gtin — (0.25 -I- z0.16) fm, we considered another parametrization listed 
in Table ITU which gives a value twice as large, namely UriN — (0.52 -I- iO.33) fm, and at the same time allows one to 
describe well the above-listed reactions. More details are given in These two sets of parameters are used to study 
the sensitivity of the result to the r]N interaction model. 

The NN interaction is included in S, P and D waves in both isospin states, namely ^5*0, "^Pj (J — 0, 1, 2), and ^£'2 
waves for T = 1 and '^S'l, ^Pi, and ^Dj {J = 1, 2, 3) for T — 0. The corresponding off-shell amplitudes were obtained 
using the separable representation of the Paris potential from |2l| . The remaining higher partial waves of the NN 
subsystem were taken as plane waves. 

The Faddeev equations for the r]NN system were solved only for the states of lowest angular momentum where a 
relative two-particle S state, either (NN) or {7]N), is coupled with an S state of the third particle relative to this 
pair, resulting in the partitions rj + {NN) or + (rjN), respectively. As the explicit calculation shows, in the energy 
region of current interest, Q < 60 MeV, these S waves are almost totally overlapping. The other partial waves were 
treated perturbatively up to first order in the 77 A^ and NN scattering amplitudes. 

The separable representation of the two-body t matrices allows one to reduce the original Faddeev equations to a 
set of two integral equations in one variable only, presented schematically in Fig. ^ Their structure is identical to 
the usual Lippmann-Schwinger equation for the two coupled channels (A^ + A^*) and (77 + d), where the notation "d" 
is used for the interacting NN subsystem. The reaction amplitude M is then expressed in terms of the auxiliary 
quasi-two-body amplitudes and Xn- as is presented in the same figure. With respect to the application of three- 
body scattering theory to reactions with three free particles in the final state we refer to • Some technical points 
concerning the r]NN problem are described in 

We begin the discussion of our results with the total cross section shown in Figs. |2| and |2j for pp rjpp and 
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np — > rjnp, respectively. As already noted, the three-body result is normalized to the data at one point, namely at 
g = 36 MeV to CT = 5.6 ^ih using PROMISE- WASA data ^ for pp and to cr = 35 fih with CELSIUS data JlJ for np. 
Of course, the ratio between the different curves is fixed by the model. Therefore, apart from the fact that the same 
normalization factor is used for both reactions, the only nontrivial result in Figs. [21 and |31 is the energy dependence 
of the cross sections. For obvious reasons, detailed experimental data are available for the pp —^ rjpp reaction only, 
which, however, is much more difficult for a theoretical study of the near threshold region because of the Coulomb 
interaction between the protons. In the present work we have not included the Coulomb force, which fact, of course, 
hampers to some extent a quantitative comparison with the data. Nevertheless, the results shown in Figs. |21 and |3| 
allows us to draw the following conclusions: 

The three-body model reproduces quite well the energy dependence of the experimental cross section in the pp 
channel. The deviation very close to threshold must primarily be ascribed to the neglect of Coulomb repulsion. As 
for the np channel, the agreement with the data of 18] appears to be less satisfactory. The theory exhibits a more 
flat form of the cross section. Clearly, for a meaningful conclusion about the quality of the present model a greater 
amount of data is needed, especially in the vicinity of the threshold. 

The first-order calculation, where the FSI is reduced to only rjN and NN rescatterings, does not reproduce the 
form of the pp — > rjpp cross section as predicted by the three-body treatment and also observed in the experiment. 
As was noted in the introduction, it is expected that the leading terms associated with rjN and NN rescatterings do 
not constitute a reliable approximation of the whole multiple scattering series, due to the slow convergence of this 
series in the vicinity of the virtual rjNN pole. The strong difference between the first order and the exact three-body 
calculation even at Q = 60 MeV demonstrates the crucial importance of higher order terms. We consider this result 
as a direct evidence for the fact that a three-body approach is definitely needed for a correct treatment of the rjNN 
interaction. Thus, the first-order approximation, although seemingly reasonable in view of the weak r]N interaction 
as compared to the NN one, is in general not relevant. 

The question of a correct treatment of FSI in rj production is closely related to the relative fraction of S waves in the 
final state. As was mentioned above, in accordance with the short range hypothesis for the 77 production mechanism, 
the S waves in the rjNN system are expected to provide the major part of the cross section in the threshold region, 
thus making the FSI effect to be of particular importance in this reaction. In this context it is instructive to compare 
the relative contribution of the S waves to the total cross sections of — > rjnp and pp — + rjpp. The results of the lA for 
both processes are presented in Fig. 21 where the ^Sqs contribution 24] is shown separately. Because of the relatively 
long range nature of the mechanism for 77 photoproduction on the deuteron, the contribution of higher partial waves 
in the final state increases rapidly and already at the energy Q = 60 MeV the S wave provides only a 10 % fraction 
of the total cross section. In contrast to this, the ^Sqs configuration strongly dominates the reaction pp rjpp in the 
whole energy region, so that even at Q = 60 MeV it produces about 77 % of the cross section. In conclusion, the role 
of the final state interaction depends not only on the phase space available for the final particles (which is the same 
for both reactions) but much more crucially on the range of the primary NN —^ NN* transition interaction. 

On the right panel of Fig. (21 the dependence of the results on the model ingredients is exhibited. Firstly, we would 
like to note a visible sensitivity of the cross section to the choice of which meson exchange governs the production 
potential. Namely, tt dominance tends to make the form of the total cross section more convex as compared to the 
full curve where p exchange prevails, whereas the lA gives in both cases practically indistinguishable results except 
for the normalization factor C„. Using a weaker rjN interaction from T^, we obtain, as expected, a weaker influence 
of FSI as is shown by the dash-dotted line in the right panel of Fig. (21 

In Fig.|31we compare our calculation for the pp invariant mass distribution with the COSY-11 and TOF ^3] data. 
At the excess energy Q = 15.5 MeV the three-body result was normalized to the experimental integrated cross section 
(7 = 3.24 ^b of Q. The distribution at Q = 41 MeV is given in arbitrary units and the theoretical full curve was 
normalized to the area enclosed by the data points. For a lower excess energy Q = 15.5 MeV one finds a reasonable 
agreement of the three-body calculation with the data, although one notes a systematic slight underestimation of 
the experimental cross section in the region above M^^ = 3.54 (GeV)^. At least some part of the deviation stems 
again from an overestimation of the pp interaction effect at the lower end of the spectrum due to the neglect of the 
Coulomb force. In view of the overall normalization, inclusion of the Coulomb repulsion would reduce the ^Sq peak 
near Mpp = 2Mn, and would enhance the other part of the distribution. 

However, at the higher excess energy Q = Al MeV the shape of the theoretical distribution is in strong disagreement 
with the data. In particular, the broad maximum around M^^ = 3.62 GeV^ is not explained. Of course, in order 
to reach at a definite conclusion, the Coulomb interaction must be included in addition, but it is very unlikely that 
this would have such a strong effect as to explain the shape. Moreover, our results are at variance with those of ^2 
where an appreciable contribution of higher angular momentum states in the pp subsystem is predicted, whereas our 
calculation docs not indicate a sizeable influence from P and D partial NN waves. The corresponding contributions 
can be estimated by comparison of the full and the dash-dotted curves in Fig. O In the latter calculation only the 
lowest ^Sqs configuration in the rjNN final state was included. 
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On the other hand, it is worth to underHne the importance of the NN tensor force, which increases quite appreciably 
the transition to the ^Si state from an intermediate ^Di after 77 production. The relevance of this intermediate ^Di 
state is related to the strong momentum transfer accompanying ij production and emphasizing the high momentum 
components of the NN wave function. One of the consequences of this fact is that the enhancement by the NN 
interaction for the T ~ transition is even slightly stronger than that for T = 1. This result is at variance with 
where a ratio of the enhancement factors of T = 1 to T = of 1.85 at relative NN momenta pr > 0.5 fm^^ was 
found. According to our results in Fig.El this ratio is reduced from 1.35 to only 0.98 at pr = 0.5 fm^^ if the transition 
^Di — >^5i is included. 

In Fig. 13 the calculated distribution with respect to the 77 kinetic energy in pp ijpp is plotted together with 
the CELSIUS data • One readily notes quite a good agreement except for a very small region close to the maximal 
kinetic energy. We do not consider this discrepancy as a serious defect because, as was noted in some fraction 
of the events must have been lost. Furthermore, we would like to point out that the data of the pp mass distribution 
at Q = 41 MeV in the right panel of Fig. appears to be in conflict with the data at Q = 37 MeV in the right panel 
of Fig. 13 The latter increases monotonically with meson energy except very close to the maximal kinetic energy T,;, 
whereas the TOF results for da/dM^p in Fig. [3 exhibit a pronounced broad maximum at intermediate pp masses. In 
view of the simple relation 

dcr da 

one would expect on the other hand a similar behaviour for both distributions, because the small difference in the 
energies can hardly result in a very different form of the spectrum. Thus new experimental data for the T,, and/or 
distributions in the region Q = 37 — 41 MeV are needed in order to resolve this discrepancy. 
Finally, we present in Fig.|S|the angular distribution of an emitted proton for two excess energies. The dominance 
of p exchange in the production mechanism leads to a pronounced isotropy of the cross sections. This becomes 
apparent by comparison with the angular distributions using a pure tt exchange. They exhibit a much stronger 
angular dependence with a more pronounced peaking at forward and backward directions, which is consistent with a 
more periferal character of the tt exchange mechanism. In general, the proton angular distribution with a weak cos^ 9p 
dependence is consistent with the almost isotropic form of the experimental results of 13J, although the observed 
magnitude at Q = 41 MeV is overestimated by about 20 %. Furthermore, one should keep in mind the pronounced 
disagreement between theory and the TOF data for the M^^ spectrum at the same energy (see Fig.lSJ. With respect 
to the rj angular distribution, we did not analyze it here in view of its obvious sensitivity to higher tjN resonances, 
predominantly £'i3(1520), which were omitted in the present calculation. 



IV. CONCLUSIONS 



In the present paper we have investigated rj production in NN collisions in the near threshold region. The main 
attention was focussed on the understanding of the role of FSI, which is expected to strongly influence the reaction 
dynamics near threshold. We believe that in conjunction with previous work P, S ^ i5j the present study was 
necessary in order to establish a unified picture of this reaction. In particular, it can be quite useful for the planning 
and analysis of future experiments. From our results we draw the following conclusions: 

(i) The three-body treatment of the rjNN interaction provides a satisfactory, quantitative explanation of the energy 
dependence of the total cross section for pp — > rjpp. However, the dependence of the results on the details of the 
NN NN* transition introduces some uncertainty. In particular, assuming a pion exchange dominance in the 
production potential leads to stronger energy dependence of the cross section and requires, therefore, a weaker -qN 
interaction in order to bring the theory into agreement with the data. 

(ii) There is a discrepancy between our results and the data of |0 for the M^^ distribution at a cm. excess energy 
Q = 41 MeV. Namely, the broad peak of the data is not described by the theory. Comparing to other theoretical 
work, we do not reproduce the result of 15| , where this broad structure is explained by the contribution of higher 
partial waves in the NN subsystem. On the contrary, according to our flndings the contribution of P and D waves 
in the final pp subsystem plays a minor role only, at least in the region up to Q = 60 MeV considered in the present 
work. 

(iii) Even close to threshold, where the final NN system is predominantly in an S wave state, the tensor NN force in 
the r = channel must be taken into account. As we have shown, the intermediate transition to a ^Di configuration 
with subsequent rescattering into the ^Si state enhances sizeably the NN interaction effect in the triplet state. 

Before concluding, we would like to note that the sensitivity of the FSI effect to the NN NN* transition 
mentioned in (i) indicates that in principle the short range part of the final state wave function must be calculated 
more thoroughly than is done here. Obviously, the most straightforward way would be a three-body treatment of 
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the whole NN — > rjNN process without restriction to the final state. However it seems to be very difficult in view 
of the great amount of channels which have to be included in such an approach, not to mention the pure theoretical 
ambiguity in treating the NN system within the three-body formalism. On the other hand, one can expect that short 
range mechanisms may be included perturbatively. For example, the t:NN configurations could be considered in first 
order only (apart from the driving meson exchange potential) by including the additional mechanism where a pion, 
created in an NN collision, produces in turn an rj via rescattering on one of the nucleons. Clearly, such an approach, 
which would allow us to adopt the same FSI concept, is practically much easier to manage than a rigorous solution 
of the three-body coupled channel model. 
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Appendix A 



In this appendix we give the explicit expressions for the amplitudes Mj"\ a £ {7r,77, p} for the transition NN — > 
NN* determining the transition matrix element TNN-iPTP')- The latter was defined in the cm. system with p and 
p' denoting the momenta of the initial and final states, respectively. In the following expressions we have introduced 
the notations: 



e = E + MN, e' = E' + MN, e* ^ E* + Mn- , 



P^ = iE + E', p + p') , P; = {E + E*, -p~ p') , AM* = Mn' - Mm , ^^^^ 
ff'i*) — £ eL f (*) — E -I- -P— 

The expressions for p exchange are obtained by using the on-shell Gordon decomposition of the nucleon current 

u{p') 7p u{p) = u{p') [(p^ +p^) - icr^^b,. - p'^)] u{p) (A2) 

and the analogous relation for the pNN* coupling. In the present paper the invariant energy lun- of the TV* or rjN 
subsystem was determined according to the spectator on-shell scheme 

u;^, = W-MM-^i^ + —^ ), (A3) 

2 Mjv Mjy + nijj 

which is consistent with the nonrelativistic three-body treatment of the rjNN system. Furthermore, this choice of 
LON* justifies the use of the Gordon decomposition for the tensor part in the pNN coupling since in this case the 
nucleon remains on-shell after emission of a p meson. The only approximation, whose validity is not justified is the 
Gordon relation for the pNN* coupling. We believe, however, that, since the coupling constant g'^ is treated here as 
a free parameter, it is a reasonable approximation. 

Using the vertices ((H)) one finds the following contributions from tt, r/, and p meson exchange: 
(i) TT meson exchange: 



-2M^y AM J, Mm, ^-^^ ml ' ^^^^ 

where 

F3(4) =0, (A6) 

^6(5) = (A7) 

^7(8) = ±-^7)- (A8) 
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(ii) 77 meson exchange: Since the formal expressions are the same as for tt exchange except for a 
couphng constants and meson mass, we do not need to repeat them here: 



(iii) p meson exchange: 



where 



M 



(p) 



2Mn V ^MnMn' - ml 



1 ' 



-^2(1) 
F3 



T(.9y +5t) 



q' ■ p W 



9v + 9t 



e 

9T 



-2- 



p ■ r 



P-P 



f' ■ P 



ee' / AAf * AAf * 



F4 = 

+ 

Fa 



2MNe 
gv +9T 



W 



e* 

9T 



ee 

'4 



p ■ r 



1 



AM 



2MNe 
9v + ffT 



e V ee 
W + e + AM* - 



J n 



AM* AM* 



p2 pMp* 



AM* 



ee' 



2e' 



5T 



+ 



2M7vee' 

9v + 9T 
ee' 

9T 



P-P 



(AM* 



2e 



.(H^ + e' + M+^( 
V e / e* V 



+ e' + AM* + — 
e 



p -q 



AM* 



2(VF + M 



A' 



AM 

„2 



„/2 



^(M^ + e + AM* + q+^(w^ + e' + ^)-p'.r^ 
e* V e / e V e* / AM* 



5v + 9t 



9v + gr 
ee'e* 



2(H^ + Mjv) 



e* AM 
2{W + Mat) 



5t 



2A/Are2e' 



l¥ + e' 



/ 2 pfj.p*. 



AM* 



AM* 



2MAree'e* 



2 ppp* 

+ e + AM* + — + — — ^ 
e AM* 



replacement of 
(A9) 

(AlO) 

(All) 

(A12) 

(A13) 

(A14) 

(A15) 
(A16) 
(A17) 
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TABLE I: Coupling constants and form factor parameters of the aNN vertices of ©. 



2 


An [GeV] 


2 

47r 


A^ [GcV] 


2ii 

47r 


2 
9t 
47r 


Ap [GeV] 


UA 


1.7 


0.4 


1.5 


0.84 


31.25 


1.4 



TABLE IL The same as in Table^for the aNN* vertices. 



9nNN* 

V3 


a; [GeV] 


g-nNN' 


a; [GeV] 


5t 


1.45 


0.404 


2.00 


0.695 


1.7 
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M 



11 

N 



d 



X 



ci d 



X,, 



N* 



N + 



N 



d 



Xn, 
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N 



Xn.> 



N* 



N 



/7I, p, T] + 

N 



Xn. 



N 
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'1 


+ 


X, 


7" 



N* 



N 



FIG. 1: Diagrammatic representation of the three-body integral equations for the reaction A''A^ r/NN. 
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FIG. 2: Left panel: Total cross section for pp rjpp. Notation of curves: dotted: impulse approximation; dashed: rescattering 
in two-body subsystems only; dash-dotted: rescattering in A'^A'' subsystem only ; full: three-body calculation. Data from 
compilation j^. Open triangle: normalization point {Q,ct) — (36,5.6) from |23l . Inset shows near-threshold region. Right 
panel: Dependence of three-body result on model ingredients. Notation of curves: solid identical to the one in the left panel; 
dashed: without p exchange in production potential; dash-dotted: calculation with Sii(1535) parametrization from TF] with 
weaker t^A'^ interaction, a,,jv = (0.25 -|-i0.16) fm, and tt dominated A'^A*' NN* potential. 
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FIG. 3: Same as the left panel of Fig.|5|but for np rjnp. Data from |^ 
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FIG. 4: Total cross sections for pp -qpp (left panel) and yd 
Dashed lines include only ^5*05 wave in the final state. 
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■qnp (right panel) calculated in the impulse approximation. 
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FIG. 5: Distribution of squared pp invariant mass in pp r]pp. Notation of curves as in the left panel of Fig. |5| The 
dash-dotted curve is obtained with a pure Sos final state. Data at Q = 15.5 MeV are from i and at Q = 41 MeV from 0. 




FIG. 6: Enhancement factors associated with the A'^A'^ interaction in T = and T = 1 channels as function of the squared 
relative A'^A'' momentum. The dashed and dash-dotted curves are obtained with and without inclusion of the ^Di — » ^Si 
transition, respectively. 
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FIG. 7: r; meson spectrum for pp rjpp as a function of rj kinetic energy T,,. Solid and dotted curves show respectively the 
three-body model and lA predictions. Data from 19]. 
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FIG. 8: Angular distribution of protons in pp rjpp. In the dashed curves only vr and rj meson exchange is included in the 
production potential. Data from |l^. 



